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Low-noise optical lattices for ultracold ®Li 
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We demonstrate stable, long-term trapping of fermionic ®Li atoms in an optical lattice with MHz 
trap frequencies for use in a quantum gas microscope. Adiabatic release from the optical lattice in the 
object plane of a high-numerical-aperture imaging system allows a measurement of the population 
distribution among the lowest three bands in both radial directions with atom numbers as low as 
7 X 10^. We measure exponential ground band heating rates as low as 0.014(1) s“^ corresponding to 
a radial ground state 1/e lifetime of 71(5) s, fundamentally limited by scattering of lattice photons. 
For all lattice depths above 2 recoil, we find radial ground state lifetimes > 1.6 x 10® recoil times. 

PACS numbers: 37.10.Jk, 42.60.By, 07.60.Pb, 42.60.Mi 


Trapped quantum particles are used widely in modern 
atomic physics from quantum information science Q, [H 
and quantum simulations of many-body physics Si to 
atomic clocks i and studies of fundamental physics Si 
i. All of these experiments benefit from long mo¬ 
tional coherence times, often because they enable coher¬ 
ent rather than statistical averaging of results. Such long 
times require preparation of the particles in a well-defined 
motional state of the trap, ideally by ground state cool¬ 
ing. The traps must not only be stable enough to prevent 
the particles from escaping, but they should preserve the 
carefully prepared state of motion for as long as possible. 

Its light mass, m, makes fermionic ®Li particular^ 
suited to quantum simulations in optical lattices [8|. 
All energy scales in an optical lattice are naturally 
parametrized by the lattice recoil energy, /iiZreci and recoil 
frequency, v^ec = h /associated with the geomet¬ 
ric lattice spacing, a, where h is Planck’s constant. For 
the same tunneling rate in recoil units, the absolute tun¬ 
neling rate is a factor of 14.5 (6.7) faster for ®Li than for 
®^Rb ("^'^K) atoms. Assuming typical atomic lifetimes of 
one minute, it will thus be possible to study thermaliza- 
tion processes and superexchange dynamics on timescales 
much longer than previously accessible Q. Here, we 
demonstrate an intensity-stable, high-power optical lat¬ 
tice for ®Li atoms. The optical lattice is designed for 
a quantum gas microscope where individual sites of the 
optical lattice and individual atoms can be resolved in 
fluorescence microscopy OiEi- 

Fluorescence imaging of ®Li with resonant light at 
Xp = 671 nm is hampered by the large resonant recoil 
energy Ep = /(2mXp) = h x 74 kHz. Each scattering 

event adds > 2E„ on average, regardless of the atom’s 
motional state |12| ■ This recoil heating makes it challeng¬ 
ing to keep the atoms cold enough to suppress tunneling 
while scattering 0(10^) resonant photons to form an im¬ 
age. For this reason, we have to combine a laser cooling 
scheme with deep optical lattices and MHz trap frequen¬ 
cies. Trap frequencies in the MHz regime require a trap¬ 
ping laser with low intensity noise because parametric 
heating rates due to laser intensity fluctuations increase 
quadratically with trap frequency Oilil. Trap quality 
can be degraded further through thermal lensing effects 
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FIG. 1. (color online), (a) Simplified ®Li level structure (not 
to scale). We use an incoherent mixture of the two lowest 
states in the ®Li ^Si /2 ground state manifold, (b) Front view 
of the glass cell. The 1064 nm optical lattices reflect off of a 
superpolished substrate at shallow angle of incidence (70°). 
Fluorescence induced by the probe beam is collected through 
a high NA imaging system whose object plane is ~10 fim 
below the substrate, (c) Top view of the glass cell and def¬ 
inition of the lab frame {X,Y,Z). The optical lattices are 
retroreflected along X and Y. 


in the lattice optics at high intensities. To address these 
challenges, we implemented a high-power and low-noise 
lattice laser system based on Yb-doped fiber amplifiers 
seeded with an intensity-stable Nd:YAG laser. 

In this paper, we use this laser system to demonstrate 
stable trapping of the two-dimensional ground band in 
the object plane of our quantum gas microscope with 
1/e lifetimes exceeding one minute. The corresponding 
heating rates are measured with a sensitive band map¬ 
ping technique The high numerical aperture (NA 

= 0.87) of our imaging system allows measurement of 
the band populations for total atom numbers as low as 
7 X 10^. We find that the measured heating rates are 
consistent with a rate equation model based on measured 
intensity-noise spectra and spontaneous scattering rates, 
which dominate the heating at high trap depths. 

The lowest two hyperfine manifolds in ®Li are shown 
in Fig. [Ija), and the states are commonly labeled |1)- 
|6) according to their energy splitting in a magnetic bias 
field. In our experiment, we load an incoherent mixture 
of atoms in states |1) and |2) into a high-power opti¬ 
cal dipole trap (1064 nm, 300 W) and transfer its fo¬ 
cus into the center of the fused silica vacuum cell shown 
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in Fig. [n^b) and (c). The atoms are then transferred 
into a crossed optical dipole trap formed by incoherent 
light derived from a superluminescent diode at 780 nm, 
whose short coherence length avoids fringing when pass¬ 
ing through the imaging optics M- The crossed dipole 
trap is located 80 pm below the object plane of a high- 
resolution microscope system and the sample is evapo¬ 
rated further in this trap at a magnetic held of 300 Gauss. 
We then load a single layer of an optical accordion lat¬ 
tice and decrease the accordion incidence angle from 
88° to 70°. The angular change simultaneously com¬ 
presses the sample and transports it to a distance of 
10 pm from the superpolished mirror that is the hnal 
lens of the imaging system. We then adiabatically load 
~3.5 X 10^ atoms into a three-dimensional optical lattice 
along directions X, Y , and Z. As shown in Fig. [IJc), 
the A = 1064 nm optical lattice beams reflect off of the 
superpolished mirror resulting in a standing wave along 
the axial direction with spacing az = 2 cose ~ 
for an incidence angle 9 = 70°. By retrorehecting each 
lattice beam, we obtain a non-interfering optical lattice 
along the X and Y axes with equal spacings ax = dy = 
= 569 nm. The radial (axial) lattice spacing corre¬ 
sponds to a recoil frequency v^ec = 25.9 kHz (3.4 kHz). 

To image the atoms in situ, we apply the probe beam 
on the D 2 transition - containing two frequencies reso¬ 
nant with both hyperhne manifolds in the ground state 
- as indicated in Fig. [T](b). The fluorescence is collected 
on an intensihed CCD camera to produce the image in 
Fig. [21(a). Note that we have increased the field of view 
of our 0.87 NA infinite conjugate ratio imaging system 
at the expense of resolution by demagnifying the image. 

We load the optical lattice by ramping up the lattice 
powers Px and Py to 0.8 W adiabatically, as shown in 
Fig mb). Each lattice beam’s power is controlled by two 
independent servos, and control over the lattice power 
is handed over automatically depending on the setpoint. 
For setpoints below Pex = 0.95 W, we use a low-power 
high-bandwidth servo. For setpoints above Pexj we use a 
high-power low-bandwidth servo. After loading the lat¬ 
tice, Px and Py are changed adiabatically to a holding 
power and the atoms are held in the corresponding lattice 
for a variable time, thoid, during which they are heated 
by lattice intensity fluctuations and spontaneous scatter¬ 
ing of lattice photons. At the end of the experiment, we 
release the atoms from the lattice using a ramp that is 
adiabatic with respect to the band gaps of the lattice, 
but fast compared to the residual harmonic conhnement 
timescales. The high-bandwidth servo allows releasing 
the atoms within 200 ps, using the ramp shown in the 
inset of Fig. mb). At the end of the ramp, the atoms are 
allowed to expand ballistically for 1.7 ms, after which we 
apply a short probe pulse to obtain the image in Fig. me). 

This image shows clear rectangular features due to the 
band edges of the radial lattice that are convolved with 
the in-trap density distribution. We model the band 
map distribution by flat rectangular features with widths 
set by the radial lattice spacings, convolved with a two- 




FIG. 2. (color online). Starting with 3.5 x 10® atoms in 
the [X, Y) plane shown in the in-trap fluorescence image (a, 
dark image subtracted), we apply the adiabatic ramp shown 
in panel (b) and obtain the band-mapping image (c). These 
images are fit (d) with a convolution of in-trap distribution 
and band map (the deconvolved band map is shown in the 
inset). The fit residuals are shown in panel (e). All images 
represent the same region, and images (c)-(e) share the same 
color bar. We show cuts through the center of each image in 
the margins. 


dimensional Gaussian distribution. The resulting fit and 
its residuals are shown in panels (d) and (e), and cuts 
through the center of each image are shown in the mar¬ 
gins. We extract the radial band populations Cij from 
the fit amplitude in the corresponding Brillouin zone, 
shown in the inset of panel (d). 

By varying the hold time in the lattice and observ¬ 
ing how the radial ground state population Coo varies 
with time, we can measure an exponential heating rate. 
The resulting heating rates for different lattice powers are 
shown in Fig. |3l We hnd exceptionally low heating rates 
and correspondingly long lifetimes of the radial ground 
band of up to 71(5) s. Even at the deepest trap depths 
and MHz trap frequencies required to implement single¬ 
site resolved imaging in a quantum gas microscope, the 
radial ground band still has a lifetime of 20 s. For the 
deepest trap depths the heating process is dominated by 
exponential heating due to spontaneous scattering of lat¬ 
tice photons (solid line in Fig. [3]). 

Because the axial optical lattice is formed by two lat¬ 
tice beams, the position of the atoms is sensitive to in¬ 
tensity fluctuations on either lattice beam. The spacing 
of the axial lattice formed by each beam depends on the 
angle of incidence 9 as argued above. If the angles of in¬ 
cidence for X and Y are mismatched by A9 = \9x — 9y\, 
uncorrelated intensity noise on Px and Py will lead to 
fluctuations of the axial trap minimum, causing fast heat¬ 
ing that has the same dependence on motional quantum 
number as spontaneous photon scattering. By includ¬ 
ing all of these processes in a rate equation model based 
on measured intensity-noise spectra and a geometric es¬ 
timate of A9 < 1.2°, we obtain the shaded region in 
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FIG. 3. (color online). Ground band heating rates Too 
as a function of total lattice power Ptot = Pr + Py. The 
marker color indicates the fractional lattice power mismatch 
(Pr — Py)/Ptot. The heating rates are obtained by fitting ex¬ 
ponential decays to the ground band populations Coo from 
band mapping images after holding atoms in the lattice for a 
variable time. The error bars indicate the statistical uncer¬ 
tainty from the fits. We calibrate the lattice trap frequency Vx 
(uy) via lattice modulation spectroscopy for different lattice 
powers Px (Py), leading to the calibrated scales shown at the 
top for Px = Py The shaded areas indicate the results from 
a rate equation model based on the measured intensity-noise 
spectra for different angular mismatch AS. The dashed line 
indicates the simulated heating rate without intensity con¬ 
trol. The solid line shows the asymptotic contribution from 
spontaneous scattering of lattice light oc y/Ptot ■ 


Fig. [21 whose lower boundary corresponds to A9 = 0. 
In the model, the radial ground state population Coo is 
calculated by integrating over the axial state populations 
of the three-dimensional trap (see Appendix |B]) . 

The model also assumes a trap depth along the radial 
and axial axes which we calculate from measured trap 
frequencies under the assumption of sinusoidal modula¬ 
tion. The trap frequencies are calibrated via lattice mod¬ 
ulation spectroscopy [l^. The modulation spectra also 
show cross-modulation peaks at the base frequency and 
confirm the effect. For lattice powers corresponding to 
trap depths below 2hvyec) the model deviates from the 
experimental data because the harmonic oscillator ap¬ 
proximation inherent in the rate equation model does 
not describe shallow lattices well. 

The lattice beams are derived from Yb-doped fiber am¬ 
plifiers (Nufern), seeded by an intensity-stable Nd:YAG 
non-planar ring oscillator laser (Innolight Mephisto). Be¬ 
cause of the high optical power requirements, the beam 
path schematically shown in Fig. |4l(a) is designed to be 
mechanically and thermally stable. At high intensities, 
many materials show thermal lensing effects that result 
in beam pointing and focusing changes on fast timescales. 
For these reasons, the beam paths (as well as the vacuum 
chamber) use fused silica optics which are less susceptible 
to thermal lensing because of the material’s high thermal 
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FIG. 4. (color online), (a) Lattice optics setup including fiber 
amplifier, collimation optics, optical isolators, and power ac¬ 
tuators as described in the main text, (b) Relative intensity- 
noise (RIN) spectra for the X lattice laser (the Y laser has 
similar features) and the seed laser taken as described in Ap- 
pendixO (c) For low lattice powers, the high-bandwidth ana¬ 
log servo suppresses intensity noise below 10 kHz. The servo 
is conservatively tuned to work across 1.5 orders of magnitude 
in setpoint. The servo gain and bandwidth decrease with de¬ 
creasing setpoint, leading to a larger contribution of noise 
from the servo electronics. The increased noise at frequencies 
corresponding to the band gaps of the optical lattice leads 
to the higher heating rates for low lattice powers in Fig. O 
At the same time, the servo reduces the noise below the seed 
laser RIN in the frequency range corresponding to the Hub¬ 
bard model parameters |rrll . 


conductivity, low thermal expansion coefficient, and low 
index of refraction sensitivity to temperature. 

Thermal effects in Tb 3 Ga 50 i 2 (TGG) optical isolators 
resnlt in a loss of isolation at high optical powers (l8l | 
which can lead to damage to the fiber amplihers from 
backreflections. To provide sufficient isolation, we use 
two stages of TGG Faraday rotators (FR) with fused sil¬ 
ica Brewster polarizers (BP) and A/2 waveplates (HWP). 
To suppress thermal variation in the isolators, the full op¬ 
tical power is always incident on the isolators. 

The optical power in the lattice beams is controlled in 
two stages. As a first stage, we use a low-bandwidth, 
high-dynamic-range actuator built from a z-cut birefrin- 
gent quartz plate attached to a galvo motor. The quartz 
plate acts as a variable waveplate (Berek compensator) 
and in combination with a polarizer becomes a variable 
attenuator with a dynamic range of ^20 dB over ^4.5° 
of rotation. 

For the second stage of intensity stabilization, we use a 
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large-active-area acousto-optical modulator (AOM, Crys¬ 
tal Technology 3080-197) as the actuator. Thermal ef¬ 
fects in the Te02 AOM crystal can also result in point¬ 
ing drifts on slow time scales. Since the optical lattices 
are focused to 80 ^m waists at the position of the atoms, 
such pointing noise would cause position fluctuations of 
the trap center. To ameliorate these effects, the lattice 
beams are focused through the AOM crystal by adjust¬ 
ing the distance between the fiber tip and the air-spaced 
fused silica collimator (Optosigma 027-0510). 

We detect the optical power on two independent, shot- 
noise-limited transimpedance amplifiers using 1064 nm 
enhanced InGaAs photodiodes (PD, Hamamatsu G8370- 
01) with gain just small enough to cover bandwidths up 
to 700 kHz. The detected photovoltages are then used as 
the input to two independent servo loops. 

For optical powers above Pex = 0.95 W, we do not re¬ 
quire fast control over the lattice depth but require that 
the laser intensity is as passively stable as possible. In 
this regime, we feed back on the angular rotation of the 
quartz plate using a slow digital feedback loop. The low- 
bandwidth (small-signal 3 dB point ^2 kHz) actuator 
ensures that we cannot write noise onto the laser at fre¬ 
quencies comparable to the trap frequencies. 

In Fig. HKb), we show laser intensity power spectral 
densities for the X lattice and seed lasers. The spec¬ 
tra are normalized to the DC power and are commonly 
referred to as relative intensity-noise (RIN) spectra (see 
Appendix . The prominent relaxation oscillation peak 
in the seed laser spectrum becomes well-suppressed when 
engaging the built-in intensity-noise eater. After seeding 
the fiber amplifiers, the intensity stability is degraded by 
10 — 20 dB for Fourier frequencies below 1 MHz. Noise 
spikes from the switching power supplies driving the fiber 
amplifier pump diodes can be strongly suppressed by low- 
pass filtering the power supplies (MPE DS26387). Acous¬ 
tic pointing noise on the fiber tip from power-supply fans 
is suppressed by removing the power supplies from the 
ampliher enclosure. 

For optical powers below Pex, we use analog feedback 
on the rf amplitude driving the AOM with an rf mixer 
and limit the servo’s bandwidth by inserting a steep low- 
pass filter at 500 kHz (LPF-B0R5). The bandwidth lim¬ 
itation ensures that no electronic noise gets written onto 
the lattice amplitude at frequencies in the MHz regime. 
For small power setpoints, the AOM servo has a nonlin¬ 
ear transfer function due to the use of a frequency mixer. 
This nonlinearity results in setpoint-dependent gain and 
reduced bandwidth for small setpoints, which is partially 
compensated by a Schottky-diode-based linearization cir¬ 
cuit in the controller. By conservatively tuning the loop, 
we are able to control the lattice depth over 2.5 orders 
of magnitude at the expense of slightly degraded noise 
spectra and lifetime for small lattice depths, as seen in 
Fig. [3l Here, such a large dynamic range is useful to ob¬ 
tain the adiabatic band mapping images. Additionally, 
the low-power servo suppresses noise at frequencies rele¬ 
vant for quantum simulation of Hubbard models El- If 


required, the servo tuning can be reoptimized for even 
longer lifetimes for lattice depths of interest. 

In conclusion, we have demonstrated stable trapping 
of ®Li atoms in the radial ground band of 1064 nm op¬ 
tical lattices spanning 2.5 orders of magnitude in trap 
depth. We measure ground band populations with sen¬ 
sitivity down to 7 X 10^ atoms and the 1/e lifetime r 
in the ground band can exceed one minute for deep lat¬ 
tices and is longer than 10 s (27rz/i.ecT > 1.6 x 10®) for 
all lattice depths above 2hi/^ec- These heating rates are 
one to two orders of magnitude smaller than in ion traps 
with comparable trap frequencies El , with smaller dis¬ 
tances to the nearest surface (10 fim here). In a three- 
dimensional optical lattice with MHz trap frequencies, 
we demonstrate radial ground band lifetimes that are 
comparable with the longest trap lifetimes measured in 
optical dipole traps and one-dimensional optical lat¬ 
tices [2l|. The heating rates are well-explained by a rate 
equation model based on the measured intensity-noise 
spectra. These spectra can be further tailored by servo 
design for application in quantum simulation experiments 
with ®Li. Optical traps with MHz frequencies enable ion- 
trap-like spectral addressability [J, are compatible 
with proximity to surfaces, and may have applications in 
achieving stron g co upling to high-quality-factor mechan¬ 
ical resonators |^ . 
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A.M., M.F.P., and C.S.C. are supported by NSF GRFP, 
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Appendix A: RIN Measurements 

To measure the RIN, we used a shot-noise limited 
transimpedance amplifier design based on an In¬ 
GaAs photodiode (Hamamatsu G8370-01), a fast opera¬ 
tional amplifier with 800 MHz gain-bandwidth product 
(OPA843), and a transimpedance gain of 330 D. To re¬ 
duce thermal drift from interference between the pho¬ 
todiode and its cover glass, we removed the photodiode 
window. Residual light contamination was attenuated by 
placing an interference filter (Semrock FF01-1020/LP-25, 
ND = 5 for visible light) before the photodiode. 

We measured RIN power spectral densities (PSDs) by 
putting the output of such a photodetector onto a Fourier 
transform precision voltmeter (SR760, varying resolution 
band width [RBW] 65 Hz below 12.5 kHz, and 500 Hz 
above) below 100 kHz or a battery powered RF spec¬ 
trum analyzer (Anritsu MS2721A, RBW = 10 Hz) above 
100 kHz Fourier frequency. The noise spectra were nor¬ 
malized to a frequency bin width of 1 Hz (assuming white 
noise in each bin), and the SR760 spectrum was con¬ 
verted from dBVrms/A/Hz to dBm/Hz, to be comparable 
with the spectrum analyzer. Such spectra were then fur- 
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ther normalized to the optical carrier power using the DC 
voltage measured with an RMS voltmeter. For RIN mea¬ 
surements, we typically apply 12 mW of optical power to 
the photodiode (resulting in 2.0 V DC signal) to reduce 
the shot noise level below the noise floor Pfloor of the RF 
spectrum analyzer (typically Pfioor ~ —155 dBm at RBW 
= 10 Hz). Battery powered operation reduces noise from 
external sources and makes the measurements compati¬ 
ble with the SR760 results, allowing us to combine data 
from both devices in the plots shown in Fig. |4l 

Appendix B: Heating Rate Model 

From the measured RIN spectra, we calculate single¬ 
particle heating rates in deep optical lattices due to 
trap intensity fluctuations and cross-modulation for mis¬ 
matched trap centers 0, [13, [IS HI- We estimate heat¬ 
ing rates from spontaneous scattering of lattice photons 
for optical lattices in the Lamb-Dicke regime from stan¬ 
dard expressions for laser cooling [iLflAl^l^. 

We then combine all heating rates in a three- 
dimensional rate equation for the probabilities to oc¬ 
cupy the motional state n = (n^, Uy^riz) of the form 

P^{t) = Y}R T14—n+An Rr^+Ar,^r^]Pr^{t)■ (Bl) 

An 

To compare our measured band mapping data against 
the heating rate coefficients, we numerically solve 
Eqn. (|B1I) . We sum over the vertical direction to 


get horizontal band populations Cij which we can then 
directly compare to the band map fit coefficients. 

Most of the heating rate coefficients depend on the RIN 
PSD. For the low power servo settings, we linearly inter¬ 
polate (in linear units) between measured spectra such 
as the ones in Fig. |4l(c). For the high power servo set¬ 
tings, we use the power-independent RIN spectra from 
Fig. |D(b). The PSDs are linearly interpolated (in linear 
units) at the Fourier frequency of interest. For this inter¬ 
polation, the Fourier frequencies of interest are multiples 
of the trap frequencies, The trap frequen¬ 

cies are calibrated against the measured optical power 
via lattice modulation spectroscopy. 

The initial distribution among trap states is assumed 
to be Boltzmann with a temperature adjusted to match 
the measured band map pictures at short times. We limit 
the state space to states with energies below the lowest 
modulation depth (along the y direction, estimated from 
Vy) and assume that an atom is completely lost once it 
is heated to states with higher quantum numbers. 

The three-dimensional rate equation is propagated 
from the initial condition and the vertically averaged 
populations Cij are fit with exponential loss curves. The 
resulting exponential decay rates are then compared to 
the experimental data in Fig. [31 

The parameter with the largest uncertainty is the mis¬ 
match between the angles of incidence of the two lattice 
beams A0 responsible for the cross-modulation heating 
contribution. With the simulation, we generate heating 
rates for several values below a conservative upper limit 
Ad < 1.2°, leading to the shaded areas in Fig. O 
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LATTICE HEATING MECHANISMS 

We present a short summary of results derived in 
Refs. (Sll-I^ and extend them to derive an expression 
for cross-modulation of overlapped but non-interfering 
optical lattices with incommensurable periodicity in the 
tight-binding limit. To ensure we account for all numeri¬ 
cal factors, we present all equations with frequencies v in 
units of Hz in the sense of cos(27ri/t). Rates and rate coef¬ 
ficients F still carry factors of 27r in the sense of exp(—Ft) 
and are presented in units of 


1. Noise-induced heating in harmonic traps 

Intensity fluctuations of the dipole trapping laser will 
result in fluctuations of the trap frequency and produce 
parametric transitions between states of the same parity. 
In the limit of small fluctuations in a tight-binding optical 
lattice we can approximate the potential for an atom with 
mass M around a lattice minimum xq as harmonic with 
trap frequency z/j. We model the intensity fluctuations 
with a small fluctuation e and write 

Vix) = Y{2Trvtf{l + e){x - Xq)^ (SI) 

In this limit, we can also confine ourselves to transitions 
between harmonic trap states n and n ± 2. The rate co¬ 
efficients for these transitions can be calculated 

in time-dependent perturbation theory as |S2l | 

TT^ 

Rn± 2 ^n = + 1 ± l)(n ± 1), 

= Rf{iyt){n -I- 1 ± l)(n ± 1). 

where S'e(2z/t) is the relative power spectral density of the 
laser intensity in (linear) units of evaluated 

at twice the trap frequency. For an infinitely deep trap, 
these rate coefficients result in an exponential increase in 
mean energy with the rate 

F, = {E)/{E) = (S3) 

If the harmonic trap center xq is subject to fluctuations 
6, we write 


V{x) = —{27T,,tnx-{xo + S)]^, (S4) 


which induces transitions between trap states of opposite 
parity. The largest rate coefficients are given by 

R,i±l^n = 4:T^‘^^i'fSs{i2t){n 

= 'K'^v'^Ss/a{vt){np]^±]^), (S5) 

= R5{vt){n+]^ ± i), 

to show that we can normalize the position noise by the 
harmonic oscillator length = hj{M2'KVt) to get a more 
intuitive result. Assuming an infinitely deep trap again, 
we find a linear increase in mean energy with rate 

{E)/{hvt) = (hi/t)-i X ATT^My^Ss{vt) 


2. Cross-modulation 


Consider the case of two harmonic traps with different 
trap centers and different trap frequencies. The total 
potential is then again harmonic with 


L(a:) =— [(27rz^i)^(a; - xi)'^ + {2T:v2f{x 

= Y(27>-) (z/i +V 2 ) 

X [(a; - xof + r(l - r){xl + x^)] 



(S7) 


where we have defined the trap frequency asymmetry 
parameter r = v\l(y\ + 1 ^ 2 ) and the new trap center 
Xq = rxi -f (1 — r)x 2 - If we allow intensity fluctuations 
Ci in each trap frequency as in Eqn. (EU, we see that the 
trap center Xg also fluctuates and find (to first order in 
Ci) relative trap frequency and trap center fluctuations 
with 


e = rei -I- (1 - r)e 2 , 

S = r(l - r)(ei - e 2 )(xi - 2 : 2 ). 


(S8) 


Clearly, the trap center will only fluctuate if both 
traps have non-zero frequency, if their centers are not 
matched and only if the intensity fluctuations are uncor¬ 
related. Assuming independent noise processes we find 
power spectral densities 

S,(2^t) = r^S,A2i^t) + (1 - .gg, 

Ss(l't) = [?'(1 - rjAxAlSeAl't) + Se^At)}, 
with trap mismatch Aa; = Xi — X 2 - 


3. Application to surface-reflected lattices 


We are now ready to apply the heating rates derived in 
the previous Sections to the special case of our surface- 
reflected optical lattices. We model the site-local poten¬ 
tial for a deep optical lattice (in the tight-binding regime) 
as a three-dimensional harmonic oscillator. Experimen¬ 
tally, we can directly measure determine the ground 
to first excited band energy difference using parametric 
heating measurements as described in the main text. The 
intensity-noise induced heating rates along the horizontal 
lattice directions x and y are 


(SIO) 


Rli^y) = 


where we have neglected the small contribution of the x 
(y) lattice beam to the other trap frequency Vy {vx)- 
Along the vertical direction, the surface lattice beams 
also form an optical lattice with spacing given by the 
laser wavelength A and the angles 6^ and Oy with respect 
to the surface. The position of the n-th lattice minimum 
away from the surface is given by Zn{0i) = A/(2sin0i). 
If we assume a small angle mismatch 9^ = 0 — Ad/ 2 
and 9y = 9 + A0/2, we find that the minima produced 
by the lattice beams differ by Az = Zn{9x) — Zn(/9y) ~ 
Zn{9)A9/ tan0. Using the arguments presented in Sec. [21 
we can immediately see that in addition to the intensity 
noise heating, the position of the n-th vertical lattice min¬ 
imum will shake as well. We hnd heating rates for the 
vertical direction as 

Rei^^) =y(2^,) + (1 - 


Rsi^z) 




(Sll) 


x[s/{u.) + sy{v^)]. 


From Eqn. (|S11|) . we can immediately see that the cross¬ 
modulation of the vertical lattice center can be problem¬ 
atic because it scales with the cube of the vertical trap 
frequency and is proportional to the intensity noise power 
spectral density evaluated at the fundamental of the trap 
frequency. However, careful angle-matching of the inci¬ 
dent lattice beams should be able to reduce this heating 
rate dramatically. 


4. Single-photon scattering 

The optical dipole trap will also produce heating from 
lattice photon scattering. Each scattering event will pro¬ 
duce an average total energy increase of rh. Erpr, w here 
r/j = 2 for an optically thin atomic cloud |S4l-IS6l|. In 
a motional state resolved picture in the Lamb-Dicke 
regime, the scattering events will produce heating that 
is described well by transitions between states n ± 1 •<— 
n to first order in the Lamb-Dicke parameter squared 


rf' = Erpr / (hvf ). This treatment leads to transition rate 
coefficients [s3, [Hi 


Rn±i^n ='^ Y??^[rsc(A) -I- rsc(A =F 27ri/t)] 
. 1 , 1 


(S12) 






with identical dependence on n as the rate coefficient 
for trap center fluctuations. Here, the scattering rate 
rsc(A) is to be understood as the steady-state scattering 
rate solution from the optical Bloch equations for detun¬ 
ing A from atomic resonance. For the far-detuned case, 
rsc(A =F 27rzzt) ~ rsc(A) = Tsc. Since we do not distin¬ 
guish between scattering events due to individual lasers, 
we added a factor of 1/3 in Eqn. (IS12|) . which assumes 
isotropic heating and ensures that the total energy in¬ 
creases as [H | 


{E) — 'y^{Ej) — hl/jR/^Pnj (t) — rhErec^sc: 

j j n=0 

. 

where we used that the populations are normalized to 
=1. 

To relate the scattering rate Tgc to experimentally ac¬ 
cessible observables, we note that the scattering rate is 
proportional to the overall optical dipole trap depth U. 
The atomic polarizability a for an alkali atom ground 
state is well described by a single Lorentz oscillator 
model |S9l | with resonant frequency ujq and linewidth T: 


a = 


r/u 


u!q — uj^ — /loq)T 


(S14) 


The proportionality coefficient between Tgc and U is 
given by the ratio of the real and imaginary parts of 
a [S9l |. and we hnd for ®Li ground state atoms in a 
A = 27rc/w = 1064 nm optical dipole trap that 





2 Ima 
Re a 


~ 1.09 X 10“®. 


(S15) 


In summary, the heating rate for lattice photon scat¬ 
tering in the Lamb-Dicke regime is 


R^ 

J. 


Ej-pc 

3 Ti huj 


(S16) 


Our best estimate for the total dipole trap depth U 
is the lattice modulation depth along the vertical direc¬ 
tion which for a sinusoidal modulation is related to the 
measured vertical trap frequency by 



(S17) 


where = h/i^SMdl) is the geometric recoil frequency 
related to the vertical lattice spacing = A/(2sin0). 
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